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Abstract 



The Coulomb corrections to the hehcity amphtudes of high-energy photon spht- 
ting are examined. The consideration is based on the amphtudes obtained exactly in 
the parameter Za within the quasiclassical approach valid for small angles between 
all photon momenta. We consider the case when the transverse momenta of both 
final photons are much larger than the electron mass. It is shown that at Za ~ 1 the 
Coulomb corrections essentially change the result for the cross section as compared 
to the Born approximation. The effect of screening is also taken into account. 



1 Introduction 



One of the most interesting nonlinear QED processes at high energy is sphtting of one 
photon into two in electric fields of atoms. The total cross section of this process does not 
decrease with increasing photon energy. First results on the observation of high-energy 
photon splitting on atoms have been obtained recently in the Budker Institute of Nuclear 
Physics [0. Theoretically this process has been investigated in 0, H, ||, || only in the 
lowest order in Za (Born approximation), Z\e\ is the nucleus charge, a = e'^/An = 1/137 
is the fine-structure constant, e is the electron charge. Though the expressions obtained 
in 0, ^ are rather cumbersome, some numerical calculations based on the results of these 
papers have been carried out in |^ P] . Using the Weizsacker- Williams method providing 
the logarithmic accuracy the cross section of the process has been obtained in an essentially 
simpler form in [Q. The comparison of the exact cross section with the approximate 
result 1^ has shown that the accuracy is better than 20%. The Coulomb corrections to the 
cross section, which can essentially modify the result as compared to that obtained in the 
Born approximation, have been unknown up to now. 

Recently we derived in the analytical expressions for the high-energy photon-splitting 
amplitudes. The result was obtained exactly in the parameter Za at small angles /2 and /s 
between the momenta k2 , ks of the final photons and the momentum ki of the initial one. 
This region of angles gives the main contribution to the total cross section of the process. 
Small angles and high energies of photons allow one to use the quasiclassical approach 
developed in p, Q at the investigation of coherent photon scattering in a Coulomb field 
(Delbriick scattering). This approach gives the transparent picture of the phenomenon and 
essentially simplifies the calculation. 

In the present paper we start from the analitycal results obtained in and investigate 
the role of the Coulomb corrections in the photon splitting process. We restrict ourselves 
to the case |k2_L| = ^^2/2 ^ \^3±\ =1^3/3 rri ( tUj = |kj|, m is the electron mass) 
when the amplitudes can be essentially simplified. A pure Coulomb potential as well as 
the influence of screening is considered. 



2 Kinematics of the process 

Below we use the coordinate system with z-axis directed along ki so that = aki/a;i 
and a_L = a — az^i/uji for an arbitrary vector a. According to the uncertainty relation 
the lifetime of the virtual electron-positron pair is r ~ |r2 — ri| ~ uji/{rn? + A^), where 
A = max(|k2_L|, Iks^l) ^ uoi. The characteristic transverse distance between the virtual 
particles can be estimated as (m^ -|- A^)~^/^, which is much smaller than the length of 
the electron-positron loop. The characteristic impact parameter is q ^ 1/A , where A = 
k2 + ka — ki is the momentum transfer. At small k2_L and ksx (/2,3 <^ 1) we have 
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The characteristic angular momentum is / ~ c<j/A ^ 1 , and the quasiclassical approxima- 
tion can be apphed. 

Let us discuss a screened Coulomb potential. In the Thomas-Fermi model the screening 
radius is Vc ~ (ma)~^Z~^^^ . If i? -C 1/A -C Tc (i? is the nucleus radius), then the 
screening is inessential and the amplitude coincides with that in the pure Coulomb field. 
At 1/A ~ Tc the screening should be taken into account. Obviously, the impact parameters 
g ^ Tc do not contribute to the total cross section. Due to this fact we shall concentrate 
ourselves on the momentum transfer region corresponding to the impact parameter g < Vc- 
If 1^2x^/^72 + ksx^/cus -C , then it follows from ([l|) that the condition Q < holds 
only when |A_|_| = |k2_L + k3_|_| > r~^. Thus, the main contribution to the amplitude is 
given by the region of momentum transfer Aj^, restricted from below. In addition, at 
Lojivp? + A^) ^ Tc the angles between ki 2,3 and ri 2,3 are either small or close to vr , and 
corresponding expansions are used in our calculations. 

According to the Furry theorem the photon-splitting amplitude is an odd function with 
respect to the parameter Za. In the Born approximation the amplitude is proportional 
to the Fourier transform of the Coulomb potential (~ Za//S?). Therefore, the region of 
very small momentum transfers A ~ is essential, and screening should be taken into 
account. In next orders of perturbation theory with respect to the parameter Za ( Coulomb 
corrections ) the integral over all momenta corresponding to the external field should be 
taken provided that their sum is equal to A. Therefore, even at A ~ r^^ each momentum 
is not small and the screening can be neglected. In the Born approximation the screening 
can be taken into account by multiplying the amplitude by the factor [1 — F(A^)], where 
F(A^) is the atomic electron form factor. Thus, to find the photon-splitting amplitude in 
a screened Coulomb field it is sufficient to solve the problem in a pure Coulomb field. 



3 Amplitudes 

It is convenient to perform the calculations in terms of the helicity amplitudes 

iA2A3(ki, k2, k3). The longitudinal components of the polarization vectors Cj can be elim- 
inated owing to the relation ejkj = which leads to = — e^kx/cij. After that within 
a small- angle approximation one can neglect the difference between vectors (62,3 )_l and 
the polarization vectors of photons propagating along ki and having the same helicities. 
Therefore, the amplitudes Mx^\2\-j,{}^i-i^2-i^z) are expressed in terms of the polarization 
vectors e and e* corresponding to positive and negative helicities, respectively. It is suffi- 
cient to calculate three amplitudes, for instance, M_| (ki,k2,k3) , M4.++(ki, k2, k3) and 

M_|_+_(ki, k2, k3) . The rest amplitudes can be obtained by substitutions: 



M+_+(ki, k2, k3) = M++_(ki, k3, k2) , M_A2A3(ki, k2, k3) = M+A^Aalkl, k2, k3) (e ^ e*) , 

where A denotes the helicity opposite to A Let us introduce two-dimensional vectors 
f2 = k2_L/co'2, f3 = k3_|_/ci;3 (|f2,3| <^ 1) and f23 = f2 — f3. We represent the amplitudes 
obtained in for |k2_L| ^ m, |k3x| ^ m in the following form 
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where q is a two-dimensional vector lying in the plane perpendicular to ki, A denotes Aj^. 
In eq. the following notation is used: 

1^2 = ^2- = u^ + e , 

U2i^3a? - t^3/t2b2 .„ ^ u;2«;3C^ - ujieh"^ uiea? + t^3/t2c2 
Vi = tO , V2 = , V3 = , (3) 

Uie UJ3K2 UJ2l^3 

a = q-A + 2K2f2 , b = q+A-2/t3f3 , 
c = q+A-2£:f23 , c = q-A + 2ef23 ■ 



Note that vectors e and e* appeared in denominators in owing to the use of the 
relation 2(ea)(e*a) = a^. 

As it was mentioned above, the photon-splitting amplitude obtained in the Born ap- 
proximation 0, ^ for arbitrary energies and momentum transfers is rather cumbersome. It 
is interesting to obtain the Born amplitude from (^. In this approximation 
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(e*,q + A) (e*,q-A) 
It is convenient to rewrite the quantities in (|^) as 



1^2-1^3 4:UJ2^^3K2K3 2 



V,= ici+'^^A - — -- f^3-^0 , (4) 

\ UJi / LUi 

T> f '^3 + e \2 AuJiiU2K3e 2 

^^2 = q A 2 — ^2 , 

\ U3 J ui 

^3= q+ A + fg . 

\ U)2 J Co-I 

After that we shift the variable of integration q in each term so that the quantities 2^i_3 
become independent of the angle of the vector q. For instance, in the terms containing Vi 
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we make the substitution q — >■ q — After passing to the variable z = exp(i0) we 

can easily take the corresponding contour integral. Taking also the integrals with respect 
to |q| and e, we get for the Born amplitudes 

2tZae\i2 X fg). 
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It follows from (|^) that ln(— ai) should be interpreted as ln(— ai + iO) = In(ai) + ivr. 
Besides, 

Li2(l + ai) = 112(1 + ai - zO) = — ln(l + ai)[ln(ai) + m] - Li2(-ai) 

o 

The result (J^) is obtained for |A_l| ^ \^z\- One can show that it remains valid in the case 
|Aj_| ~ IA2I if the expression (|l]) for A^ is used in (H). Actually, in eq. (|]) the difference 
between A^ and A^ is essential only in the overall factor 1/A^. For a screened Coulomb 
potential the amplitudes (|) should be multiplied by the atomic form factor (1 — F{/S?)). 
For the case of Moliere potential |]10| it reads 
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where 

ai = 0.1 , a2 = 0.55 , = 0.35 , Pi = Pobi , (7) 



b^ = Q , 62 = 1-2 , h = 0.3 , /3o = mZ^/yi21 . 
Remind that the representation (j^) is vahd when |k2_L| , \^3±\ ^ m. 

Let us consider the asymptotics of the amphtudes (|]) at |A_l| ^ \p\, where = 
u;2f2 + uJs^s and p = {(^2^2 ~ ^3^3) f'^- It is this region of variables which gives the main 
contribution to the cross section in the Weizsacker-WiUiams approximation. To get this 
asymptotics we multiply T in (|^) by 

i = %o-q') + ^(q'-%') , 

where |A| ^ go ^ Then, for the term in @ proportional to ^{Qq — q^) one can 
put q = and A = in T and integrate by parts over q. After that, using the relation 
Vq'(9(go— q^) = ~2q5(gQ— q^) one can easily take the integral over q since at |q| = go ^ I A| 
one has 

V|q-A|; q2 

As a result, in the region |q| < go the term proportional to Za is independent of go and 
the terms of next orders in Za are small in the parameter |A|/go. 

For the term proportional to '(9(q^ — q^) we get 



|q — Z\| / |q| 



We put A = in T and perform the integration first over the angles of q and then over |q|. 
As a result, the main in go/|p| contribution is independent of go and proportional to Za. 
Taking the sum of the contributions from these two regions and performing the integration 
over the energy e, we get 
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In the small-angle approximation (|f2|, Ifs] <^ 1) the cross section of the process reads: 



^2 

da = —^\M\^ x{l- x)dxdi2d^3 , (9) 

where x = 0^2/(^1, so that = uJi{l — x). In terms of the variables p and A the cross 
section has the form 

,,^,0 dAdpdx , , 
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Substituting @ into (10) and performing the elementary integration over the angles of 
vectors A and p, we come to the expression 

AZ^a^ dp^ dA^ dx 



da 
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where the function g{x) for different polarizations has the form 
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Formulae (0) and (|T2D are in agreement with the corresponding results of |@, obtained in 
the Weizsacker- Williams approximation. However, this approach does not allow to obtain 
the amplitudes (H) themselves. The large logarithm appears after the integration of (0) 
over A^ from A^-^ up to where A^m ~ for the screened Coulomb potential and 
Amin ~ for the pure Coulomb case. It is interesting to compare the contributions of 

different helicity amplitudes to the cross section at A ^ 0. In Fig. 1 the function g{x) is 
shown for different helicities as well as the quantity 



9{x) = 9+—ix) + ^+++(x) + 9++-ix) + 9+-+ix) 



(13) 



which corresponds to the summation over the final photon polarizations. It is seen that 
g{x) has a wide plateau. 

The Coulomb corrections to the photon-splitting amplitude at A ^ are small com- 
pared to the Born term (|]). We consider the asymptotics of the Coulomb corrections at 
A in the next Section. 



4 Coulomb corrections 



To analyze the Coulomb corrections we make the further transformation of the expression 
(I). Let us multiply right side of (H) by 

1- Jdy6(y-^^)={c,^ + A^)jfi-6{{ci-A/yr-A\W-l)) , 

\ q / _^ \y\ 

and make the substitution q ^ q + A/y. After that the integral over |q| becomes trivial, 
and the integral over the angle of q can be easily taken by means of the residue technique. 
Finally, we get 



6 



M 



Ae^Za f dy signy 



Re 



1+y 



iZa 



(e*f3)(e*f: 



23, 



0^2 

j deK2 (e*, K2f2 - A) 



i-yj 



C(;2 ^ ^^3 
V f 2 ^ f 3 , 



(14) 



8(e*f3)(/s:2 + /s:2 



UJ2 \ \/s'^ — Aq'^ii^ 



(eu) s — y/s'^ — Aq'^u^ — 4(e*r)(eu) 



+ 



[Kl + 4)K2{l + l/y){e*i2){eA) 



(ef23)[£(e*r)(ef23) - ^2(1 + l/i/)(eA)(e*f2)] 



+ 



UJ2 



de 



111^ + 



UJ2 ^ UJ3 
^ f 2 ^ f 3 ^ 

2 I 



8(ef23)(/€i + £ 



Aq'^ui-s'i J \ie*ui) + i^Aq'^ul - si- 4(er)(e*ui) 

i4 + e')K2il + l/y)iei2){e*A) 
(e*f3)[K3(er)(e*f3) - /€2(1 + l/y)(e*A)(ef2)] 



+ 





UJ2 J de K3 

-^3 



feA) 



.^i(e*A) \ Mq^ul-sj 



8(e*f23)/s:3(/s:2 + £2 



^^'^iJ si - iJ4g2uf - sf - 4(e*ri)(eui)_ 



+ 



feAl 



S2 



W3(e*A) y^s2_4^2„2 

,2/, ^2 I ^2 



II + 



8(e*f2)/«3(/€i + e' 



'S2 - 4g2u2 
'd{rf — q 



2 2\1 \'^2 + ^ 



leu2j S2 + - 4g2u2 - 4(e*ri)(eu2) 
1 

+ 



+ 



2^2(e*ri) V«^2(e*ri)(ef2) - n^il/y - l)(eA)(e*f3) 

£(e*ri)(ef23) - K,il/y - l)(eA)(e*f3) , 
Here we use the following notation: 



u = A - - 1 + 

\y ^2 , 



U2 = A 1 
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Since the function R in (|I^) is independent of the parameter Za the Coulomb corrections 
can be obtained from (|14D by the substitution 
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The asymptotics of M^^^ at A — * depends on the photon helicities. The most simple way 
to get this asymptotics is to start directly from @. For M+l_ the main contribution is 
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determined by the region q p ^ A (remind that p = {uj2i2 — ^^3^3)/'^), while for M^l_^_ 
and M^l_ it comes from the region q ^ A <^ p. After the corresponding expansion and 
integration over the energy we get at A ^ 
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It follows from (16) that in this limiting case the Coulomb correction m'^^__ is small, 
while M_^4.+ and M\j^_ depend only on the direction of vector A, but not on its module 
(it becomes obvious after the substitution q qA). We discuss the role of Coulomb 
corrections in the next Section. 



5 Cross section 



As it was suggested in 1 11 1 , to overcome the problems of background in the measurement of 
photon splitting one has to register the events with |f2,3| > /o where /o -C 1 is determined 
by the experimental conditions. Let us consider the cross section integrated over fs for 
Ifsl > /q. It is interesting to compare the exact ( in Za ) cross section da/dxdi2 with that 
obtained in the Born approximation {das/dx di2) and also with the cross section in the 
Weizsacker- Williams approximation {daw/dx di2). Remind that daw/dx di2 is in fact the 
Born cross section calculated within logarithmic accuracy. Large logarithm corresponds to 
the contribution of the region A <C p = \uj2i2 — ^-^3f3)/2|, where /s ~ 3^/2/(1 — x). Taking 

up to Agjj, where (see [§) 
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and summing over the final photon polarizations we get for a pure Coulomb potential 



da 



w 



dx (if 2 



SZ^a^ g{x) ,^ f2{l-x 



In 



/2 



X 



1 — X 



/2 - /o) • 



(17) 



For the case of a screened Coulomb potential the approximate cross section is 
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The function 7 in eq. (|^) is 

7 = 1 - E «,'(lna, + 1) - , a, = 6,^ + A^J/^o' (19) 

and the coefficients a^, 6j and /5o are defined in eq. (0). If A^j„//5o ^ 1 then 7 = 
- ln(A^Ln//?o) and eq. (HD turns to eq. (0). If ^l,J 01 « 1 then 7 = -0.158. 



For the case of a pure Coulomb potential the dependence of ctq "^do jdxdi^ on /2//0 is 
shown in Figs. 2-4 for /o = 10~^, Z = 92 and different values of x, 
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In these figures the solid curves represent the exact cross sections, the dashed ones are 
Born results [daBj dx di'i) and the dash-dotted are obtained in the Weizsacker- Williams 
approximation. At x = 0.7 (Fig. 2) the difference between dasj dx di2 and da^jdxdi^ 
is small. At x = 0.3 (Fig. 3) daB/dxdi2 differs noticeably from daw/dx di2. Note that 
within a good accuracy the cross section das/dx at x = 0.3 agrees with that obtained from 



eq. (p!7|) after the integration over f2. It should be so, since da /dx is invariant with respect 
to the substitution x — >• 1 — x and at x = 0.7, as we pointed out above, the approximate 
result (p!?!) is in accordance with the exact one. At x = 0.5 a big difference between 
das/ dx di2 and daw / dx df2 (see Fig. 4) in the region f2 ~ /o can be explained as follows. 
The large logarithm appears as a result of integration with respect to fs over the range 
1(1 — x)f3 + xf2| -C x/2. After the integration over the azimuth angle ip between vectors 
f2 and — fs we should integrate over fs from /o up to x/2/(l — x) and from x/2/(l — x) 
to infinity. If x/2/(l — x) ~ /o then the contribution of the first region vanishes and 
the cross section becomes approximately two times smaller (in accordance with Fig. 4). 
In all cases the exact cross section at Za ~ 1 is noticeably smaller than the Born one. 
The magnitude of this effect depends on kinematics. For instance, at x = 0.3 (Fig. 3) 
and /2 < /o(l — x)/x when A ~ p the exact cross section is several times smaller than 
the Born one, while for /2 > /o(l — x)/x the difference of these cross sections is about 
15%. For X = 0.5 and 0.7 (Figs. 3,4) this difference is about 20% and almost independent 
of /2//o- All indicated relations between cross sections take place also when screening 
is taken into account. Emphasize that at A ~ p and Za ~ 1 the exact cross section 
differential over all variables {da/dx di2di3) is much smaller than the Born one. In Fig. 
5 this differential cross section is shown for the case of a screened Coulomb potential. 
The peak for azimuth angle = vr corresponds to small A. There is a narrow notch at 
/s = 3^/2/(1 — x) which corresponds to the condition Aj^ = 0. The width Sf^ of the notch 
is about max{Az/uj3, Po/us). For the parameters used in Fig. 5, Sf^ is about 10""^. 

Let us discuss now the cross section da/dx. In the Weizsacker- Williams approximation 
for a pure Coulomb potential the cross section daw/dx obtained from (|l^) is 
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If lui/q/^I — x) <^ Po, then the corresponding cross section in a screened potential has 
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the form 



daw 
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X \ Pq I 
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The inequahty A <^ p which provides the apphcabihty of the Weizsacker-Wilhams ap- 
proximation corresponds to a small angle (p between the vectors £2 and — fs ( when the 
vectors £2 and fs have almost opposite directions). So, it is interesting to consider the quan- 
tity daB{(pmax)/dx which is the Born cross section integrated over the angle ip from —ipmax 
to (fmax- In the case of a pure Coulomb potential the dependence of {'k f^ao)'^ daB^Pmax) I dx 
on ipmax is shown in Fig. 6 for different x and /o = 10~^. One can see that the cross section 
is saturated at relatively large Pmax- The same conclusion is valid for the case of a screened 
Coulomb potential. 

Let us represent the exact cross section as a sum da /dx = das/dx + dac/dx. As it 
was mentioned above, the region of small momentum transfers A is not important for the 
Coulomb corrections dac/dx. Therefore, one can neglect the effects of screening in the 
calculation of dac/dx. The quantity F = f^ao)'^ dac / dx is independent of /o and uji 
and is the function of Za and x. The dependence of F on x for different Z is shown in Fig. 
7. On can see that the Coulomb corrections diminish the cross section of the process. To 
realize the magnitude of this effect we plot in Fig. 8 the cross sections da/dx, das/dx and 
daw/dx for Z = 83 and /o = 10"'^. It is seen that the Coulomb corrections are important. 
This is the consequence of the fact that at Za ~ 1 the constant added to the logarithm, 
i.e. the quantity x"^ F/g{x), is large enough in a wide range of x. Although this quantity is 
independent of /o and ui, the relative contribution of the Coulomb corrections to the exact 
cross section depends on these parameters since the large logarithm contains them. 

Due to the gauge invariance the cross section da/dx should be equal to zero at x = 1 if 
the electron mass is not neglected. This is not the case for massless particles as it has been 



noted in |T2|. That is why the cross section da/dx calculated in zero-mass limit does not 



vanish at a; — 1 (see Fig. 8) . When x ^ 1 the main contribution to the cross section comes 
from the range of angles /s ~ /2/vT^-a; and /2 ~ /o , A ^ ^^2/2, where the Weizsacker- 
Williams approximation is not applicable. Since /a ^ 1, the relation 1 — a; 3> /g should be 
fulfilled. Besides, the condition k^,^^':^ m means that (1 — x)uji ~ y/1 — xcui /o ^ m or 
1 - a; > (m/^i/o)2. 

If we represent the amplitude M of the process as a sum of Born amplitude Mb and 
Coulomb corrections Mc, then jMp = |Mijp + 2 Re(M^Mc) + |Mcp. Taking into account 
in Mc only the lowest in Za term (proportional to (Za)^), we get dac^ /{Za)"^ = ci + 
C2{Zay , where ci and C2 are independent of Za and C2 > . The dependence of 
{Za)^^{7cfQao)~^ dac/dx on Za is shown in Fig. 9 for different x. One can see that 
starting from Z ~ 15 the contribution of the next order Coulomb corrections essentially 
modifies the behavior of dac. 

Thus, the Coulomb corrections are very important for the adequate description of high- 
energy photon splitting and must be taken into account at the comparison of theory and 
experiment. 

We are grateful to V.S. Fadin and E.A. Kuraev for useful discussions. 
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Figure captions 



Fig. 1. Function g{x) from eq. (0) for different polarizations: (x) (dotted 

curve), 5f++_(x) (dashed curve), 5'+++(x) (dash-dotted curve), and g{x) (sohd curve), eq. 

Fig. 2. (7Q^da/dxdi2 versus /2//0 for a pure Coulomb potential, /o = 10~^, x = 0.7, 
Z = 92, (Jo is given in the text. The dash-dotted curve corresponds to the Weizsacker- 
Williams approximation, the dashed curve gives the Born approximation, the solid curve 
is the exact result. 

Fig. 3. Same as Fig. 4 but for x = 0.3 . 

Fig. 4. Same as Fig. 4 but for x = 0.5 . 

Fig. 5. Differential cross section da/dxdf2df^ versus /s in a screened Coulomb potential 
for different azimuth angle between vectors ¥2 and fa; Z = 83, x = 0.7, ui = IGeV, 
/2 = 5mrad. The dashed curve (Born approximation) and the solid curve (exact cross 
section) correspond to = vr. The dash-dotted curve (Born approximation) and the dotted 
curve (exact cross section) correspond to = 7r/2. 

Fig. 6. The dependence of {nfQcro)~^do'{(pmax)/dx on (pmax for a pure Coulomb potential 
at different x: x = 0.5 (dotted curve), x = 0.7 (dashed curve), and x = 0.9 (solid curve); 
/o = 10-^ 

Fig. 7. The dependence of {n fQao)'^ dac / dx on x for Z=32 (dotted curve), 47 (dash- 
dotted curve), 64 (dashed curve) and 83 (solid curve). 

Fig. 8. The dependence of {nfQaoj^^da/dx on x for Z=83. The dotted curve corre- 
sponds to the Weizsacker- Williams approximation, the dashed curve gives the Born ap- 
proximation and the solid one shows the result exact in Za. 

Fig. 9. The dependence of {7ifQao)^^da/dx/Za'^ on Z for x=0.5 (sohd curve) and 
x=0.8 (dashed curve). 
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Fig. 5 
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